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EVALUATION 


The  temporal  and  spatial  occurrence  of  auroral  clutter  have  been 
described  for  an  HF  backscatter  radar.  The  clutter  signals  have  been 
characterized  by  their  intrinsic  reflectivities  which  permit  the  evalua¬ 
tion  of  models  of  auroral  clutter  which  are  developed  independently. 

This  effort  contributes  to  Technical  Planning  Objective  1C  Surveillance 
by  providing  a  better  understanding  of  the  properties  of  the  polar  iono¬ 
sphere  and  technology  for  system  simulation  and  diagnosis  in  the  presence 
of  ionospheric  irregularities. 

The  signal  processing  techniques  developed  under  this  effort  provide 
the  necessary  methodology  and  associated  software  to  serve  as  the  basis  for 
the  RADC/EEP  support  for  the  testing  phase  of  the  Experimental  Radar  System 
that  is  currently  under  development  by  the  CONUS  OTH-B  41 4L  SPO. 

RING  / 

Project  Engineer 
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1.  INTRODUCTION 


The  objective  of  this  study  is  to  estimate  the  III'  radar  volume  back- 
scatter  reflectivities  from  li  and  1-  region  irregularities  in  the  Polar 
ionosphere.  The  primary  source  of  data  is  the  Polar  l-ox  II  experiment, 
in  which  a  high  power  HF  backscatter  radar  was  operated  in  northern  Maine 
from  November  1971  through  November  1972.  The  coverage  area  included  the 
North  magnetic  pole,  with  the  auroral  oval  passing  through  the  coverage 
region  at  a  mean  range  of  1500  km.  The  present  study  makes  use  of  data 
from  the  2200-0800  UT  time  interval  during  the  months  of  September  and 
October  of  1972. 

Ionospheric  clutter  is  a  natural  phenomenon  that  can  degrade  the 
performance  of  radar  systems  designed  for  high  latitude  coverage.  The 
sources  of  ionospheric  radar  clutter  in  the  E  and  P  regions  of  the  high 
latitude  ionosphere  are  magnetic  field  aligned  reflecting  structures  in 
the  electron  density  field.  These  structures  exhibit  motion  which  imparts 
a  Doppler  shift  to  scattered  radar  energy.  The  backscatter  cross  section 
of  the  region  of  ionospheric  irregularities  illuminated  by  an  OTH  radar  may 
be  very  large,  making  the  detection  of  targets,  against  the  resulting 
clutter  background,  difficult  or  impossible. 

The  physical  mechanisms  producing  the  irregularities  in  the  high  lati¬ 
tude  E  and  F  layers  are  thought  to  be  different  even  though  the  occurrence 
of  both  kinds  of  irregularities  can  be  correlated  with  the  same  solar  and 
geophysical  phenomena.  Irregularities  responsible  for  E  layer  auroral 
clutter,  originating  at  altitudes  of  about  110  km.,  are  commonly  understood 
to  be  generated  by  plasma  instabilities  driven  by  the  horizontal  electrojet 
currents  of  the  auroral  E  region  [Greenwald,  1973]  .  One  of  the  first 
studies  attempting  to  develop  a  quantitative  theory  of  equatorial  E  region 
irregularities  based  upon  the  two  stream  plasma  instability  is  due  to 
Farley  [1963],  where  the  energy  producing  the  irregularities  is  supplied 
by  the  equatorial  electrojet.  Further  developments  in  plasma  instability 
theory  were  stimulated  by  the  shortcomings  of  the  Farley  theory,  leading  to 
the  gradient -dri ft,  or  E  x  _B  cross  field  instability  [Rogister  and  D'Angelo, 
1970;  Farley  and  Balsley,  1973],  and  an  adaptation  of  these  equatorial  zone 
theories  to  auroral  latitudes  [Ginzburg,  1970;  Ginzburg  and  Rubadze,  1972] . 
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The-  radar  aurora  is  related  to  the  visual  aurora,  but  significant 
differences  are  observed  between  the  two  phenomena.  Both  a  discrete  and  a 
diffuse  type  of  auroral  clutter  are  observed.  The  discrete  type  occurs 
mostly  at  nighttime  and  is  characterized  by  large  velocities,  while  the  diffuse 
type  is  normally  observed  during  daytime  (afternoon)  and  exhibits  slower 
velocities.  It  seems  that  the  diffuse  radar  aurora  is  related  to  stable 
visual  auroral  arcs;  and  the  discrete  radar  aurora  is  related  to  pulsating 
rays  and  moving  patches  in  the  visual  aurora  [Bullough  and  Kaiser,  1955]. 

Several  models  of  I!  region  irregularities  arising  from  plasma  instabili¬ 
ties  have  been  invoked  to  explain  the  range  of  ti  region  clutter  phenomena, 
particularly  its  frequency  dependence.  There  has  not  been  as  much  progress 
in  explaining  high  latitude  T  layer  clutter.  One  useful  method  of  studying 
high  latitude  !•'  layer  irregularities  has  been  by  observing  scintillations 
in  electromagnetic  radiation  originating  from  radio  stars  and  satellites. 

These  scintillations  are  produced  by  electron  density  fluctuations  associated 
with  F  layer  irregularities,  and  permit  observations  to  be  made  on  the  geo¬ 
physical  distribution  of  these  irregularities  [Aarons,  1973;  Nielson,  and 
Aarons,  1974;  Martin  and  Aarons,  1977] . 

/Another  method  of  observation  is  the  use  of  II.  F.  radar  backscatter 
which  constitutes  the  subject  of  this  report .  Due  to  the  strong  magnetic 
field  alignment  of  ionospheric  plasma  irregularities,  their  radio  scatter 
cross  section  is  highly  anisotropic,  with  a  maximum  value  in  the  plane 
orthogonal  to  the  local  geomagnetic  field  direction.  In  the  case  of  a 
backscatter  radar  this  means  that  appreciable  radar  returns  are  received 
only  when  the  propagation  direction  is  closely  perpendicular  to  the  geo¬ 
magnetic  field  direction.  Most  of  the  radar  studies  of  high  latitude 
irregularities  have  used  VI I F  and  U1IF  transmissions  which  are  insensitive 
to  I  layer  field  aligned  irregularities  because  there  is  not  sufficient 
refraction  to  direct  the  rays  perpendicular  to  the  magnetic  field  at  those 
altitudes.  At  III-'  the  ionospheric  refraction  is  sufficiently  great  to  permit 
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orthogonality  over  a  large  area  of  radar  illumination.  This  advantage, 
however,  is  achieved  at  the  expense  of  frequent  uncertainty  about  the  path 
of  the  radar  signal,  with  consequent  difficulty  in  accurately  determining 
the  location  and  size  of  the  scattering  region.  This  difficulty  leads,  in 
this  report,  to  a  simplification  in  computing  the  scattering  volume,  with  a 
resulting  uncertainty  in  the  value  of  the  reflectivity. 


2.  DESCRIPTION  OF  EXPERIMENT 

The  Polar  Fox  II  experiment  included  the  operation  of  an  HF  backscatter 
radar  by  the  Raytheon  Company  in  northern  Maine  from  November  1971  through 
November  1972.  The  coverage  area  is  shown  in  Figure  1.  The  experiment  used 
a  high  power  radar  (800  kW  peak)  to  provide  coverage  over  a  90°  azimuthal 
sector  from  -30°  to  +60°  azimuth  measured  from  true  North.  The  coverage  area 
included  the  north  magnetic  pole,  the  auroral  oval  and  an  area  at  the  most 
easterly  azimuths  which  could  reasonably  be  expected  to  be  equatorward  of 
the  auroral  zone.  The  auroral  oval  passes  through  the  coverage  region  at  a 
mean  range  of  1500  km. 

The  radar  operated  at  30  pulses  per  second  with  a  10  kHz  chirp  bandwidth 
corresponding  to  15  km  range  resolution.  Radar  parameters  are  listed  in 
Table  1 . 


TABLE  1 

Radar  Parameters 


Location : 

Frequency : 

Power : 

Range  Resolution: 
Receive  Beamwidth: 
Azimuth : 
Measurement : 


Caribou,  Maine  4  7°N  68° IV 

6  to  26  MHz 
64  Kw  Avg. 

15  km 

9  degrees  at  10  Ml lz 

11  steps  from  -30°T  to  +60T 

One  minute/Frcq . /Beam  Position/llour 


3 


4 


Figure  1.  Polar  Fox  II  Radar  Coverage  Map 


The  Polar  Pox  II  radar  used  a  4  element  vert ical ly-polarized  log- 
periodic  transmitting  antenna  array  and  a  28  element  vertically  polarized 
log-periodic  receiving  area  with  a  nominal  30  dB  taper.  Antenna  pattern 
measurements  were  performed  in  a  cooperative  program  by  Raytheon,  Lincoln 
Laboratory  [Nichols,  1973]  and  the  Rome  Air  Development  Center  [Parry,  1973] 
using  a  variety  of  aircraft.  Results  of  all  the  measurements  were  analyzed 
at  Lincoln  Laboratory  to  prepare  net  gain  tables  for  each  antenna  referenced 
to  the  receiver  input  or  transmitter  output  including  all  line  and  network 
losses.  Gain  vs.  elevation  angle  tables  were  prepared  for  each  antenna 
for  each  operating  frequency,  for  each  steer  angle  and  for  each  transmitting 
system  element  configuration  used.  Measured  transmit  and  receive  patterns 
for  3  different  frequencies  on  the  antenna  bo  resight  (15°  true  North)  are 
shown  in  Figure  2. 

Radar  operations  were  scheduled  for  two  twenty-four  hour  runs  each  week. 
Three  backscatter  ionograms,  one  on  boresight  and  one  to  east  and  west  of 
bores ight  were  made  each  hour.  The  ionograms  covered  6  to  26  MHz  in  freq., 
and  1000  to  4000  km  in  range  which  is  the  full  operating  range  of  the  radar. 
The  three  ionograms  required  about  four  minutes  to  complete.  -These  were 
followed  by  a  series  of  fixed  frequency  backscatter  surroundings  to  examine 
the  signal  strength  and  Doppler  spectrum  as  a  function  of  range  and  azimuth. 
Three  fixed  frequencies  were  selected  well  in  advance  using  LSSA  prediction 
techniques,  to  provide  backscatter  coverage  of  the  1000-4000  km  range  in 
1000  km  increments.  The  selected  frequencies  could  be  changed  based  on 
operating  experience  and  using  the  backscatter  ionograms  for  guidance. 

Operations  began  on  the  lowest  frequency.  Data  were  collected  at  a 
given  azimuth  for  one  minute,  including  a  one  second  noise  measurement  with 
the  transmitter  off,  and  average  Doppler  spectra  computed  at  each  range 
interval.  At  the  end  of  each  minute  the  average  Doppler  spectrum  at  each 
range  interval  was  written  on  magnetic  tape.  I'he  spectra  consisted  of  32 
values  at  15/16  Hz  intervals.  This  measurement  was  repeated  in  10°  azimuth 
increments  through  the  90°  scan  sector.  This  eleven  minute  sequence  was 
then  repeated  using  the  middle  frequency  selected  to  provide  ground  coverage 
in  the  2000-3000  km  range,  and  finally  the  sequence  was  repeated  at  the 
highest  radar  frequency  for  coverage  in  the  3000-4000  km  range.  The  back- 
scatter  measurement  took  thirty-three  minutes  and  was  repeated  once  each 
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Figure  2.  Polar  Fox  II  Transmit  (Tx)  and  Receive  (Rx) 
Antenna  Patterns  for  3  Frequencies 


hour  resulting  in  a  one  minute  measurement  at  a  particular  frequency  and 
azimuth  each  hour. 


3.  DATA  REDUCTION 

3.1  Footprints 

The  range/Doppler  information  recorded  in  each  azimuth  measurement 
during  the  fixed  frequency  mode  of  operation  was  displayed  as  a  footprint  as 
in  Figure  3.  This  is  a  format  in  which  the  Doppler  spectrum  of  the  received 
signal  is  displayed  in  each  of  the  IS  km  range  bins.  The  footprint  is  a 
method  for  viewing  a  three-dimensional  array  in  a  two-dimensional  format, 
with  the  amplitude  of  the  received  power  in  -dBW  in  the  appropriate  range/ 
Doppler  location  as  the  third  coordinate.  In  another  footprint  display  the 
numerical  values  of  the  amplitude  are  replaced  by  letter  symbols  representing 
10  dB  increments  to  aid  in  the  visual  discrimination  of  amplitude  levels. 

These  displays  are  used  for  the  routine  mode  analysis  described  below.  A 
three-dimensional  plot  display,  shown  in  Figure  4,  requires  a  substantial 
amount  of  computer  time  to  produce,  and  is  valuable  on  a  limited  basis  when 
amplitude  variations  are  small  or  when  clutter  patterns  are  complicated. 

They  are  also  valuable  in  making  similar  clutter  patterns  among  a  number 
of  footprints  stand  out  clearly. 

Only  the  data  from  the  fixed  frequency  mode  of  operation  were  available 
for  processing.  Recordings  of  the  swept  frequency  ionograms,  produced  once 
each  hour  just  before  the  fixed  frequency  scans,  were  available  in  poor 
quality  Fax  records  and  were  of  limited  use. 

3.2  Threshold  and  Interference 

A  threshold  level  was  established  for  each  footprint  based  upon  the 
noise  spectrum  of  the  data.  A  histogram  of  the  signal  levels  in  a  footprint 
for  all  range/Doppler  cells  showed  that  the  number  of  these  noise  signals 
typically  decreased  with  amplitude  for  about  (>  dB  above  the  minimum  signal 
level.  This  was  taken  as  the  population  of  the  noise  unrelated  to  ionospheric 
or  ground  clutter,  so  that  the  threshold  was  set  b  dB  above  the  minimum 
signal  level . 
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Interference  from  other  users  of  the  HI-'  band  is  frequently  present 
and  particularly  so  in  the  nighttime  hours  when  the  useable  band  of  fre¬ 
quencies  is  contracting  and  operators  are  competing  for  clear  channels. 

Since  the  radar  was  restricted  to  specific  operating  frequencies  there  was 
a  significant  amount  of  interference  contamination  in  the  data.  This 
interference  usually  produced  energy  in  a  narrow  Doppler  band  (typically 
2  Hz)  with  constant  received  power  at  all  ranges,  and  was  often  of  signifi¬ 
cant  amplitude  relative  to  the  clutter. 

An  algorithm  was  designed  to  routinely  eliminate  these  interference 
lines  from  the  footprints,  without  raising  the  noise  threshold,  with  a  minimum 
effect  on  the  overlapping  clutter.  It  made  use  of  the  constant  amplitude 
with  range,  and  established  when  the  population  of  the  minimum  signal  in  a 
Doppler  cell,  taken  over  all  range  bins,  exceeded  a  certain  critical  level, 
and  was  also  a  critical  number  of  dB  above  the  noise  threshold.  These 
critical  parameters  were  determined  empirically.  An  example  of  the  appli¬ 
cation  of  the  interference  algorithm  is  shown  in  Figure  5.  There  were 
initially  three  interference  lines  and  it  will  be  noticed  that  while  there 
is  some  degradation  of  the  broad  clutter,  the  line  at  -1  to  0  Hz  has  been 
quite  clearly  eliminated. 

Another  algorithm  was  designed  to  recognize  ground  clutter  in  the 
footprints,  a  process  that  was  greatly  facilitated  by  the  previous  removal 
of  interference  lines  in  the  region  of  0  Hz  Doppler.  The  algorithm  makes 
a  parabolic  fit  to  five  adjacent  amplitudes  centered  around  the  peak  amplitude 
in  each  range  bin.  The  curvature  of  the  parabola  is  evaluated  and  examined 
along  with  the  Doppler  frequency  of  the  peak.  If  the  curvature  is  above 
an  empirically  determined  critical  value,  appropriate  to  the  narrow  ground 
spectrum,  and  if  the  peak  is  within  ±1  Hz,  the  signal  in  the  range  bin  is 
designated  as  ground  clutter.  When  the  ground  clutter  occurs  along  with 
spread  clutter  in  a  range  bin,  it  is  eliminated  by  substituting  the  average 
of  the  clutter  values  at  ±3  Hz  before  the  spread  clutter  is  processed 
further.  This  is  illustrated  in  Figure  6. 

3 . 3  Mode  Identificat ion 

With  a  range  coverage  of  025  to  4000  km,  the  clutter  returns  in  the 
backscatter  data  show  the  effects  of  a  variety  of  ionospheric  conditions. 
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Figure  5b,  Interference  Lines  Removed 


lustration  of  Ground  Clutter  Removal 


They  include  cases  of  direct  E  and  l;  returns,  li  and  1;  region  backscatter 
after  one  or  two  ground  reflections,  as  well  as  ground  backscatter.  In 
order  to  minimize  the  ambiguity  in  the  identification  of  ground/ionosphere 
propagation  modes,  it  was  decided  to  limit  consideration  to  returns  from 
less  than  1500  km  range.  This  would  assure  a  high  probability  of  selecting 
only  direct  returns  from  F.  or  F  region  clutter  sources. 

The  footprint  listings,  shown  in  Figure  3,  were  the  basic  data  presenta¬ 
tion  used  for  mode  identification.  These  footprint  displays  had  a  threshold 
of  minimum  signal  level  established,  and  any  interference  lines  were  re¬ 
moved.  Mode  identification  was  based  essentially  upon  a  visual  recognition 
of  more  or  less  distinct  clutter  regions  in  the  range  interval  from  625  km 
to  1500  km.  A  distinct  region  of  clutter  confined  to  less  than  1100  km, 
the  line  of  sight  horizon  range  for  an  altitude  of  110  km,  was  designated 
as  E  clutter.  A  distinct  region  of  clutter  commencing  beyond  this  range 
would  be  designated  as  F  clutter.  In  the  example  shown  in  Figure  3,  the 
close  in  returns  (less  than  820)  are  very  likely  from  the  E  region.  The 
broad  spectrum  returns  commencing  at  1225  km  and  continuing  to  nearly  2000  km 
are  designated  as  F  returns.  The  narrow  spectrum  returns  centered  at  0  Hz 
from  the  longer  ranges  are  ground  clutter.  Usually  a  region  of  clutter 
commencing  at  less  than  1100  km  and  extending  well  into  the  F  region  range 
without  apparent  discontinuity  would  all  be  designated  as  F  clutter.  These 
fairly  simple  rules  provided  the  basis  for  scanning  a  large  quantity  of 
data  to  produce  statistically  accurate  samples  of  E  and  F  clutter  returns 
from  which  to  obtain  estimates  of  reflectivities. 

Median  values  of  the  reflectivities  from  designated  E  and  F  sources 
did  indeed  seem  to  indicate  that  these  two  sets  of  reflectivities  consti¬ 
tuted  two  distinct  populations.  A  test  was  devised  to  check  whether  these 
populations  represented  two  different  clutter  sources,  or  whether  the 
distinct  populations  could  have  resulted  from  the  different  assumptions 
regarding  layer  height  with  consequent  antenna  gain  and  ionospheric  loss 
values,  made  in  the  computations.  Using  a  five  day  sample  of  data, 
reflectivities  of  returns  which  had  been  designated  as  E  layer  clutter 
were  computed  using  F  layer  assumptions  -  i.e.,  F  layer  height,  and  cor¬ 
responding  antenna  gain  and  ionospheric  absorption  loss.  The  resulting 
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population  was  clearly  different  in  median  value  from  that  based  on  the 
assumption  that  the  reflectivities  represented  and  li  layer  source.  This 
test  increased  confidence  that  the  above  mode  identification  criteria  led 
to  a  statistically  valid  discrimination  of  propagation  modes,  and  therefore 
clutter  sources. 


4.  RADAR  EQUATLON 

The  volume  backscatter  reflectivity  is  intended  to  be  independent  of 
the  specifications  of  the  measuring  system  and  to  characterize  the  intrinsic 
nature  of  the  auroral  backscatter.  For  the  volume  reflectivity,  it  is 
necessary  to  define  a  volume  in  the  ionosphere  which  contributes  to  the 
received  power  at  a  specific  time  delay.  The  azimuthal  beam  width  and 
the  pulse  width  define  two  of  the  dimensions.  The  vertical  dimension  is 
defined  as  the  angle  between  two  rays  orthogonal  to  magnetic  field  lines 
and  separated  by  a  range  difference  equal  to  the  pulse  length.  This  is 
depicted  in  Figure  7  and  was  determined  by  three  dimensional  ray  tracing 
using  model  ionospheres  to  have  a  nominal  value  of  0.5°. 

The  volume  reflectivity  is  derived  from  the  radar  equation  as 
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Illuminated  volume  = 


(R0)  (RAc)  m3 

i  =  pulse  width 

(A  =  azimuthal  beamwidth 

At:  =  effective  vertical  angle 

The  loss  term  was  adapted  from  a  method  developed  for  ionospheric 
absorption  calculations  [George  and  Bradley,  1974]. 

The  auroral  absorption  loss  term  was  determined  from  the  model 
developed  by  Foppiano  [Foppiano,  197S;  Vondrak  et  al,  1978]  of  monthly  median 
absorption.  This  is  a  model  of  the  spatial  and  temporal  variation  of  the 
vertical  aurora  absorption  mapped  in  terms  of  a  coordinate  system  of 
corrected  geomagnetic  latitude,  longitude  and  local  time.  The  locations 
of  D  region  penetration  for  H  and  F  height  reflections  were  determined  in 
the  geomagnetic  coordinate  system  and  used  to  estimate  the  vertical  absorp¬ 
tion  at  30  Ml  1  z .  Tlie  zenith  angles  of  the  rays  were  calculated  to  modify 
the  vertical  absorption  values  by  a  secant  factor  and  the  values  were 
normalized  to  the  frequency  of  interest  through  the  inverse  frequency- 
squared  dependence  of  non  deviative  absorption. 

3.  RFSULTS  OF  ANALYSIS 

The  returns  analyzed  here  were  collected  during  the  months  of  September 
and  October,  1972,  between  the  nours  of  2200  and  0800  UT.  Of  the  nine 
days  of  operation  during  the  month  of  September,  data  from  one  day  could 
not  be  used  because  the  transmitter  power  was  not  recorded.  The  system  was 
also  in  operation  nine  days  during  the  month  of  October,  and  only  the  data 
from  several  isolated  time  segments  could  not  be  processed  due  to  diffi¬ 
culties  in  reading  the  data  tapes. 

Ihe  data  from  both  months  were  analyzed  to  study  the  volume  backscatter 
reflectivity  as  a  function  of  layer  height  and  frequency.  The  October  data 
were  analyzed  in  more  detail  in  an  attempt  to  observe  any  geomagnetic  or 
auroral  effects  on  the  reflectivity.  In  particular,  an  attempt  was  made 
to  observe  any  correlation  of  reflectivity  with  geomagnetic  coordinates  or 
with  auroral  conditions  relating  to  t he  auroral  oval. 


Using  DSMP  satellite  photographs,  four  of  the  nine  (.lays  of  operation 
in  October  were  designated  "active",  while  the  remaining  five  days  were 
designated  "quiet".  On  the  active  days,  photographs  of  the  auroral  zone 
showed  auroral  boundaries  consistent  with  a  Q  between  3  and  4.  On  the 
quiet  days  the  auroral  boundary  in  the  photographs  was  consistent  with  a 
Q  between  0  and  1.  liach  of  these  data  sets  was  further  separated  into 
evening  (2200-0200  UT)  and  midnight  (0300-0800  UT)  data.  Thus,  there  were 
four  separately  processed  sets  of  data. 

Figure  8a.  shows  the  region  covered  by  the  Polar  Fox  11  radar,  with 
the  southern  auroral  oval  boundaries  for  active  and  quiet  conditions  in 
the  local  evening  time  segment  superimposed.  Under  quiet  conditions  the 
boundary  misses  the  closest  range  segment  and  divides  the  farthest  two 
segments  into  eastern  and  western  azimuthal  sectors.  The  active  auroral 
boundary  covers  a  good  portion  of  the  closest  range  segment  at  the  western 
azimuths.  Figure  8b.  shows  that  under  active  conditions  all  range  segments 
are  covered  by  the  oval  in  the  midnight  time  segment. 

Each  of  these  sets  of  data  was  analyzed  for  range  and  azimuth 

variation  which,  with  the  separate  time  segments,  constituted  the  basis 
for  an  attempted  correlation  with  geomagnetic  coordinates.  At  eacli  azimuth 
a  median  reflectivity,  with  upper  and  lower  quartiles,  was  computed  for 
each  of  the  three  range  groups  of  F  layer  returns,  and  for  the  one  range 

group  of  E  layer  returns.  The  range  analysis  was  intended  to  show  any 

strong  correlation  with  geomagnetic  latitude,  while  the  analysis  of  evening 
and  midnight  time  segments  was  intended  to  show  any  strong  correlations 
with  geomagnetic  longitude.  No  attempt  was  made  to  resolve  the  azimuthal 
variation  into  these  geomagnetic  coordinates. 

5 . 1  Geomagnet ic  Dependence  of  Reflectivity 

A  comparison  of  simultaneous  UHF  observations  of  auroral  clutter  and 
high  latitude  scintillations  at  F  region  altitudes  [Unger,  1974 ]  indicate 
that  E  region  and  F  region  clutter  sources  may  have  different  geomagnetic 
latitude  distributions.  It  is  found  that  auroral  clutter  peaks  considerabl 
south,  geomagnetical ly,  of  the  most  intense  scintillations.  On  the  other 
hand,  a  comparison  of  scintillation  observations  at  the  same  F  region 
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altitudes  with  the  optical  aurora  observed  in  I  (MSP  satellite  photographs 
[Martin,  1977]  shows  a  strong  geomagnetic  correlation.  The  IIP  backseatter 
data  reported  on  here  provide  simultaneous  observations  of  the  same  type 
of  both  E  region  and  F  region  ionospheric  clutter,  perhaps  making  it 
possible  to  better  infer  a  relationship  between  the  two  clutter  sources. 

5.1.1  E  Reflectivity 

Since  the  E  layer  returns  originated  in  only  the  first  range  interval, 
no  range  variation  in  the  E  reflectivity  is  presented.  Figure  9  shows 
the  E  reflectivities,  under  combinations  of  evening  and  midnight,  and 
active  and  quiet  conditions,  plotted  against  azimuth  at  three  frequencies. 
There  is  some  increase  in  reflectivity  at  all  three  frequencies,  and  a 
pronounced  increase  at  8  MHz,  toward  the  western  azimuths  under  quiet 
midnight  conditions.  This  may  be  consistent  with  the  fact  that  the  oval 
partially  covers  the  western  sector  of  the  closest  range  segment.  There 
may  be  a  similar  tendency  during  active  evening  periods.  There  was  no 
data  available  at  8  MHz  for  active  evenings,  nor  for  the  eastern  azimuths 
at  12  MHz.  Under  quiet  evening  and  active  midnight  conditions  the  oval 
either  completely  misses  or  completely  covers  the  closest  range  segment. 
Nevertheless  there  are  still  tendencies  for  the  reflectivity  to  increase 
toward  the  western  azimuths.  From  these  data,  then,  any  evidence  that 
the  intensity  of  the  irregularities  increases  in  the  auroral  oval  is 
ambiguous. 

It  may  be  observed  in  Figure  9,  showing  the  number  of  E  layer  returns 
plotted  against  azimuth  for  each  frequency,  that  under  active  midnight 
conditions  the  number  of  events  increases  toward  the  west  at  all  frequencies, 
while  no  clear  trend  exists  under  any  of  the  other  sets  of  conditions. 

Also  Figure  10  shows  some  consistent  tendency  for  the  E  layer  reflectivi¬ 
ties,  averaged  over  azimuth,  to  be  highest  under  active  midnight  conditions. 

5.1.2  F  Reflectivity 

In  the  F  layer  reflectivity  data  shown  in  Figures  10a.  and  10b.,  no 
consistent  azimuthal  trend  is  apparent  in  any  range  group  under  any  set 
of  conditions.  However,  there  is  an  increase  in  reflectivity  with  range, 
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under  almost  all  conditions,  with  the  largest  increase  occurring  from  the 
closest  to  the  middle  range  at  8  and  12  MHz .  This  seems  to  he  consistent 
with  the  supposition  that  any  range  increase  would  be  toward  regions  of 
more  intense  irregularities.  But  an  examination  of  the  auroral  ovals  show 
that  some  azimuthal  trend  in  the  range  increase  of  reflectivity  would  he 
expected,  especially  where  the  oval  covers  only  a  western  portion  of  the 
range  sectors,  if  there  were  a  correlation  between  the  auroral  oval  and 
the  intensity  of  the  irregularities. 

In  the  absence  of  consistent  azimuthal  trends  in  reflectivity,  con¬ 
clusions  about  the  range  dependence  of  the  reflectivity  are  based  upon 
azimuthal  averages.  It  is  clear  from  Figure  11,  which  shows  azimuthal 
averages  of  the  reflectivity  vs.  range  for  both  the  September  and  October 
data,  that  there  is  a  consistent  increase  in  reflectivity  with  range  at  all 
frequencies,  as  well  as  a  consistency  in  the  range  behavior  at  each 
frequency.  But  without  any  ordered  azimuthal  behavior  it  is  difficult  to 
make  any  correlation  with  geomagnetic  coordinates,  or  to  make  a  geophysical 
interpretation . 

Other  explanations  than  geophysical  ones  may  be  considered  for  the 
range  variation  in  the  F  region  reflectivity.  It  was  mentioned  above  that 
there  was  considerable  uncertainty  in  the  ray  geometry,  leading  to  the 
assumption  of  1 ine-of-sight  propagation  for  computing  reflectivities.  In 
particular,  any  error  in  the  scattering  volume  from  this  assumption  would 
show  up  as  an  incorrect  range  dependence  of  the  reflectivity.  But  studies 
of  the  effective  volume  using  ray  tracing  in  model  ionospheres  indicate 
that  the  scattering  volume  at  the  shortest  ranges,  particularly  just  beyond 
penetration,  are  underestimated  by  straight  line  propagation .  The  radar 
equation  shows  that  this  would  lead  to  larger  rather  than  smaller  values 
for  the  reflectivity. 

Consideration  was  also  given  to  the  possibility  that  a  poor  estimation 
of  reflection  heights,  giving  rise  to  erroneous  elevation  angles  and 
therefore  to  erroneous  antenna  gains,  could  have  a  systematic  influence 
on  range  dependence  of  the  reflectivity.  Since  very  few  elevation  angles 
would  be  above  that  of  the  peak  antenna  gain  at  8  MHz,  and  since  the 
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dominant  error  was  to  underestimate  the  reflection  height  and  thus  under¬ 
estimate  the  elevation  angle  and  antenna  gain  for  shorter  ranges,  any 
systematic  error  in  the  reflectivity  would  he  on  the  high  rather  than  the 
low  side. 

Also  it  seems  possible  that  returns  from  the  closest  ranges  could  he 
scattered  back  from  nonorthogonal  reflections.  A  consistent  off- 
orthogonality  of  the  closest  F  layer  returns  due  to  insufficient  refrac¬ 
tion  could  lead  to  a  reduced  amplitude  at  the  shortest  ranges.  Estimates 
of  the  decrease  in  amplitude  for  off -orthogonal  reflections  run  to  6  dB/deg. 
Table  2  shows  reflectivities  of  F  layer  clutter  listed  against  frequency 
and  range.  The  reflectivities  are  listed  separately  for  each  of  the  four 

sets  of  conditions  being  considered.  In  each  listing  the  change  in 

2  3 

reflectivity  with  frequency,  in  m  /m  /km  is  listed  for  each  frequency. 

This  shows  that  change  of  reflectivity  with  range  is  larger  at  the  lower 
frequencies,  contradicting  the  expected  behavior  for  the  above  orthogonality 
related  variation. 


5. 2  Frequency  Dependence  of  Reflectivity 

While  a  spatial  variation  in  reflectivity  is  usually  interpreted 
as  corresponding  to  the  spatial  distribution  of  the  irregularities  producing 
the  clutter,  a  frequency  dependence  should  correspond  to  an  intrinsic 
property  of  the  clutter  source.  The  frequency  dependence  relates  theoreti¬ 
cally  to  the  relative  intensity  of  irregularities  of  different  sizes.  A 
given  wavelength  of  incident  radiation  is  reflected  according  to  the  presence 
of  corresponding  size  irregularities.  Much  theoretical  work  has  been  done 
to  relate  the  structure  of  the  irregularities,  shown  in  the  frequency 
dependence  of  the  scattered  radiation,  to  the  plasma  instability  models 
proposed  to  explain  the  irregularities  [R.  A.  Greenwald,  1975;  M.  D'Angelo, 
1975].  However  the  models  have  been  developed  initially  to  explain 
irregularities  generated  in  the  E  region  by  the  equatorial  electrojet;  and 
attempts  have  been  made  to  extend  them  to  the  auroral  electrojet.  There 
has  been  much  less  progress  in  developing  models  of  F  region  irregularities, 
especially  at  high  latitudes.  Some  current  two-stream  and  gradient  drift 
models  predict  a  power  spectrum  for  the  scattered  radiation  of  the  form 


P'k  n,  where  k  is  the  magnitude  of  the  wave  vector,  with  n  =  3.-3. 5 
[Chestnut,  1968;  Leadabrand,  Larson  and  Hodges,  1967).  However,  some 
measurements  show  a  wider  range  of  values  for  this  power  law  exponent 
depending  on  conditions  and  the  frequency  of  the  scattered  radiation. 

5.2.1  E  Reflectivities 

For  the  purpose  of  studying  the  frequency  dependence  of  reflectivity, 
both  the  E  layer  and  F  layer  reflectivities  were  averaged  over  azimuth. 

The  frequency  dependence  of  the  E  layer  reflectivities  is  shown  in  Figure  1 
for  each  of  the  four  conditions  established  above.  While  the  active 
midnight  reflectivities  are  somewhat  higher  than  those  under  other  condi¬ 
tions,  the  trend  in  all  cases  is  very  similar.  The  power  law  exponent  for 
these  reflectivities  has  the  exceptionally  high  value  of  16,  as  opposed  to 
the  theoretical  value  on  the  order  of  3. 


Figure  12.  E  Layer  Reflectivity  vs.  Frequency 
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The  l:  layer  reflectivities  shown  in  Figure  13  for  September  and 
ictober  are  averaged  over  the  four  sets  of  conditions  as  well  as  over 
azimuth.  Here  the  l‘J  MHz  data  fjr  October  was  included  to  compare  with 
corresponding  data  from  September.  There  was  no  10  MHz  data  in  October. 

'Ihe  frequency  dependence  of  the  i  layer  reflectivity  is  very  similar  for 
the  two  months  with  no  consistent  differences. 

i  Reflectivities 

Table  1  shows,  in  addition  :o  the  range  variation  in  reflectivity 
referred  to  above,  the  variation  in  reflectivity  with  frequency.  These 
quantities,  estimates  of  the  exponent  in  the  power  law,  are  listed  under 
each  column  of  reflectivities  at  constant  range.  Again,  no  systematic 
variation  with  range  or  geomagnetic  location  is  apparent.  These  values 
are  on  the  order  of  those  computed  for  the  E  layer  clutter.  F  reflectivities 
averaged  over  range  and  the  different  sets  of  conditions  are  plotted  against 
frequency  in  Figure  14  for  the  months  of  September  and  October.  The  value 
of  the  exponent  for  September  is  again  about  16,  with  the  October  value 
being  slightly  smaller.  This  figure  also  shows  a  difference  of  6  dB  or 
less  between  the  median  l:  and  F  reflectivities  at  any  of  the  5  frequencies. 

(» .  SUMMARY 

Estimates  have  been  computed,  from  a  large  sample  of  data,  of  HF  volume 
backseat  ter  reflectivities  as  a  function  layer  height,  geomagnetic  coordinates, 
auroral  conditions  and  frequency.  F  layer  reflectivities  were  computed  as 
a  function  of  azimuth,  range  and  local  time.  While  there  was  a  strong  range 
dependence  at  all  azimuth  positions  the  range,  azimuth  and  time  variations 
c  >uld  not  be  clearly  coordinated  with  a  set  of  geomagnetic  coordinates.  The 
source  of  the  similar  range  vari;tion  at  all  azimuths,  with  no  apparent 
geomagnetic  correlation  is  not  clear.  A  number  of  possible  causes  related 
to  the  radar  system  and  analysis  procedures  were  investigated  and  were  found 
to  be  unlikely.  E  layer  reflectivities  were  computed  in  only  one  range 
interval.  These  showed  a  somewh; t  more  positive  correlation  with  geomagnetic 
coordinates  and  auroral  conditions,  tending  to  be  larger  at  those  positions 


28 


FREQUENCY  (MHz) 


covered  by  the  auroral  oval.  The  magnitudes  of  the  h  and  F  layer  reflectivi¬ 
ties,  connected  for  both  1)  region  and  auroral  absorption,  were  very  close 
for  the  same  range  interval. 


Azimuthal  averages  of  11  and  I  layer  reflectivities  were  compared  for 
the  months  of  September  and  October.  The  frequency  variation  of  both  the 
K  and  F  layer  reflectivities  showed  a  close  correspondence  for  the  two 
months,  as  did  the  range  variation  in  the  F  layer  reflect i vi t i es. 

The  analysis  presented  here  is  probably  weakened  by  the  lack  of  detailed 
ionospheric  information.  This  necessitated  the  as sump t i on  of  straight 
line  propagation  in  computing  ray  paths  and  reflection  heights  and  volumes. 
While  this  renders  the  values  of  the  reflectivities  uncertain,  no  unexpected 
systematic  behavior  in  the  reflectivities  could  be  traced  to  this  source. 
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